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ABSTRACT 

For decades, the determination of the mean density of matter (CIm) has been tied to the distribution 
of light. This has led to a "bias," perhaps as large as a factor of 2, in determining a key cosmological 
parameter. Recent measurements of the physical properties of clusters, cosmic microwave background 
(CMB) anisotropy and the power spectrum of mass inhomogeneity now allow a determination of Qm 
without "visual bias." The early data lead to a consistent picture of the matter and baryon densities, 
with il B = 0.039 ± 0.0075 and Q M = 0.33 ± 0.035. 

Subject headings: cosmology: cosmological parameters, early universe 
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1. INTRODUCTION 

The mean mass density of the Universe is a cosmolog- 
ical parameter of great importance. It can be expressed 
as a fraction of the critical density, Qm = Pm/Pci-h, or 
in physical units (gem -3 ). Moving between the Hubble 
constant: 

P m = 1.88(0 M /i 2 ) x 10" 29 g cm" 3 (1) 

where as usual, Hq = 100/ikmsec -1 Mpc -1 

gem 

Over the past thirty years much effort has been devoted 
to determining the mean matter density (see e.g., Faber & 
Gallagher, 1979; Dekel et al, 1997; Turner, 2000; Primack, 
2001). The task is daunting: the density of matter in a 
large enough sample of the Universe to be representative 
(> 10 6 Mpc 3 ) must be determined. The most well devel- 
oped techniques have been tied in one way or another to 
the distribution of light; e.g., mass-to-light ratios and pe- 
culiar velocities (see e.g., Bahcall et al, 1995; Carlberg et 
al, 1997; Dekel, 1994; Willick & Strauss, 1998). 

The mass-to-ratio technique begins from the deceptively 
simple equation that relates the mean luminosity density 
and the mean mass density: 

PM = (M/L)C. (2) 
Using the canonical value for the B-band luminosity den- 
sity today 1 and solar units, Eq. 2 can be written as the 
familiar, 

M = (M/L) B /I200h. (3) 

The task is now transformed to determining the mean 
mass-to-light ratio. Either (M/L) must be measured in a 
large enough volume to be representative (> 10 6 Mpc ), 
or measured for a system that can be argued to be "typ- 
ical." Needless to say, the accuracy of the inferred mean 
matter density can be no better than that of the mean 
mass-to-light ratio. 

In attempting to determine {M/L) attention has been 
focused on clusters, because of their size and relatively well 

1 Important questions still remain about the mean luminosity density, 



determined total masses. However, clusters are not large 
enough to provide a representative sample of the Universe, 
and because the cluster environment is a high density one 
(10 2 -10 3 times the mean density) in which only a few per- 
cent of galaxies find themselves, there is no a priori reason 
to believe that their mass-to-light ratio is representative. 

The CNOC sample provides the largest and best under- 
stood sample of clusters, and from it Carlberg et al (1997) 
infer = 0.19±0.04. By comparing cluster data with N- 
body simulations Bahcall et al (2000) find a similar value, 
Q M = 0.16 ±0.05. 

Field galaxies are more representative, but problem- 
atic because galaxy halos are large. In fact, no definitive 
evidence yet exists that the total mass of even a single 
galaxy halo has been determined. The advent of weak- 
gravitational lensing to probe the far reaches of galaxy 
halos and to determine their total masses may soon sur- 
mount this hurdle (see e.g., Fischer et al, 2000). However, 
for now, the values for VLm inferred from clusters (quoted 
above) represent the state of this art in mass-to-light ra- 
tios. 

Peculiar velocities and bulk flows are directly tied to 
the underlying mass distribution; comparison of peculiar 
velocities and mass inhomogeneity (Sp/p) fixes the mean 
density (see e.g., Dekel 1994). However, it is the distribu- 
tion of galaxies - not mass - that is measured in redshift 
surveys. For this reason, this technique probes = fl° A f /b, 
where b is the assumed linear bias factor between light 
and mass. (Nonlinear effects break the degeneracy, but 
there has been little success in exploiting this fact.) The 
peculiar-velocity technique has its own convergence prob- 
lem: Unless the mass distribution is surveyed to a volume 
sufficiently large to include all of the inhomogeneity that 
gives rise to the peculiar velocity, an accurate value for Q,m 
will not be obtained. Currently, there are unexplained dif- 
ferences in the values obtained for (3, which range from 0.4 
to almost 1. 
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A first step toward a physically based technique was pro- 
posed by White et al (1993). The method goes like this. 
Assume that galaxy clusters provide a "fair sample" of 
matter in the Universe, i.e., that their composition reflects 
the universal baryon-to-total mass ratio (= Qb/^m)- 
Then, from our knowledge of the baryon density from big- 
bang nucleosynthesis (BBN), the mean matter density can 
be inferred: 

Qm = (xt) ^(BBN) (4) 

\ M B J c i uste r 

While the bias factor b never enters, there are potential 
systematics: Can one find all the cluster baryons? Can 
the ratio be determined on large enough scales to ensure 
that baryon settling is not an issue? White et al (1993) 
and others (e.g., Mohr et al, 1998) have argued that the 
answer to these questions is yes, and several recent recent 
reviews of the mean matter density have given significant 
weight to this technique (see e.g., Dekel et al, 1997; Turner, 
2000; and Primack, 2001). 

Today, other physical measurements, including CMB 
anisotropy and measurements of the power spectrum of 
matter inhomogeneity, can be added to the mix. As I will 
describe, the physically based methods have become ma- 
ture enough to self-consistently and robustly pin down the 
mass density with a relatively small error bar. The value 
I derive from the current data is: 

n M = 0.33 ± 0.035 (la) (5) 

There are however several potential sources of systematic 
error (see below). Nonetheless, I believer that this marks 
an important first step toward an accurate determination 
of the mean matter density that is independent of the re- 
lationship of mass to light. 

2. THE INPUT DATA: PHYSICAL MEASUREMENTS 

A suite of physical measurements can now probe the 
distribution of matter and the mean matter in an unbi- 
ased way. They include accurate determinations of CMB 
anisotropy on angular scales down to 0.1 degrees and mea- 
surements of the shape of the power spectrum of matter in- 
homogeneity. CMB anisotropy depends significantly upon 
the physical baryon and matter density, as well as other 
parameters, but not upon the connection between mass 
and light. 

The shape of the power spectrum of mass inhomogene- 
ity depends significantly upon Um/i and Qm/Qb] it can 
be measured by weak-gravitational lensing (cosmic shear) 
and through the distribution of galaxies (redshift surveys). 
Cosmic shear measurements are insensitive to bias, but are 
much less mature. Redshift surveys are sensitive to bias 
in determining the power spectrum itself. However, the 
shape of the power spectrum is only sensitive to strong 
scale-dependent bias on large scales around the bend in the 
power spectrum associated with matter-radiation equality 
(~ 30Mpc/n M h 2 ). 

With MAP, Planck, full results from the SDSS and 2dF, 
and additional cosmic-shear measurements coming, the fu- 
ture for an accurate, physically based measurement of SIm 
is very bright. I believe there is now enough data to make 
a preliminary estimate of Om using these techniques. 

The following physical measurements comprise my input 
data for determining the baryon and total matter densi- 



ties: 

Power Spectrum f^M^ 
CMB Anisotropy Q, M h 2 

n B h 2 

BBN n B h 2 

Clusters Qb/SIm 
Hubble constant h 



0.20 ±0.03 
0.15 ±0.07 
0.16 ±0.04 
22+ 004 

U.UZZ_Q QQ3 

0.020 ±0.001 

(0.07±0.007)/i" 3/2 (Xray) 
(0.08±0.01)/i~ 1 (S-Z) 
0.72 ±0.07 



The shape of the power spectrum of matter inhomo- 
geneity depends upon Qm/i and Q,b/^m, as well as other 
parameters. Currently, the best determination of these 
two parameters come from an analysis of 160,000 rcdshifts 
in the 2-degree Field Galaxy Redshift Survey (2dFGRS; 
Percival et al, 2001). The baryon-to-total matter ratio is 
determined from the presence of "baryon bumps" in the 
power spectrum; it is hardly a significant result at the mo- 
ment and gives almost no weight to the determination of 
VLm and Early results from the Sloan Digital Sky 

Survey (SDSS) are consistent with the 2dFGRS numbers, 
fl M h = 0.19 ± 0.04 (Dodclson et al., 2002; Szalay et al., 

2001) . 

Significant improvement in both il,Mh and Qb/^m can 
be expected when the 2dF and SDSS surveys have amassed 
and analyzed their full samples of 250,000 and 600,000 red- 
shifts respectively. The results will also be more robust. 
For example, current analyses assume scale-invariant per- 
turbations (n = 1) and there are still correlations between 
as and Qj^h. While cosmic-shear measurements of large- 
scale structure are not yet good enough to have significant 
leverage on Q,m and fls, they will be in the future (see 
e.g., van Waerbeke et al., 2001; Maoli et al., 2001). 

The structure of the acoustic peaks in the CMB angular 
power spectrum depends upon the total matter density 
and the baryon density (in physical units): the ratio of 
the odd to even peaks pins down the baryon density and 
the height of the first peak is sensitive to the total mat- 
ter density (see e.g., Hu et al, 1997). I have used the 
values extracted from the first-year DASI results (Pryke 
et al, 2002). (The correlations between the two densities 
are small.) A re-analysis of the 1998 flight of Boomerang 
yields similar values (Netterfield et al, 2002). Combining 
the two results - which I have not done - would reduce 
the error bars by about V2. 

While the consistency of DASI, a radio interferometer 
operating at around 30 GHz, and Boomerang, bolome- 
ters in bands above 100 GHz, is very reassuring, there are 
still significant issues. Both groups used a limited num- 
ber of parameters (7) and priors on the Hubble constant 
(h > 0.45) and optical depth to last scattering (r < 0.4). 
The parameters preclude a large contribution by tensor 
(gravitational wave) perturbations; while unlikely, relax- 
ing that assumption can significantly change the inferred 
baryon density (see e.g., Wang, Tegmark and Zaldarriaga, 

2002) . With MAP and Planck on the horizon, significant 
improvement, both in the precision and the robustness of 
the results, lies ahead. 

At present, the most precise determination of the physi- 
cal baryon density comes from combining measurements of 
the primeval abundance of deuterium (see e.g., O'Meara 
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et al, 2001) with accurate theoretical predictions of the 
light-element abundances (see e.g., Buries et al, 2001). 
I follow the analysis of Buries et al (2001) in adopting 
Sl B h 2 = 0.02 ± 0.001. Here too, there are still issues to be 
resolved: possible unidentified systematics in the determi- 
nation of the primeval deuterium abundance; reliable val- 
ues for the primeval 4 He and 7 Li abundances to compare 
with the values predicted from the deuterium-determined 
baryon density. Over the next decade, as more deuterium 
systems are found and the 4 He and Li abundances are 
better understood, the accuracy and reliability of the BBN 
baryon density should improve. 

The baryon-to-matter ratio in clusters can be deter- 
mined by x-ray measurements alone and by measurements 
of the Sunyaev - Zel'dovich (SZ) distortion of the CMB 
combined with x-ray measurements. For x-ray measure- 
ments alone, I adopt the cluster baryon fraction deter- 
mined from a sample of 45 clusters by Mohr et al (1998), 
and for the SZ/X-ray determination I use the cluster sam- 
ple of Grego et al (2001). 

Potential sources of systematic error remain. There is 
some clumping of cluster gas, which could lead to an over- 
estimation of the amount of gas. The baryon-to-total mass 
ratio does vary with radius, and the calibrations are done 
by comparison with numerical simulations. The treatment 
of gas dynamics in these simulations still involves approxi- 
mations and assumptions. Better simulations, more obser- 
vations of clusters with the SZ technique and x-ray should 
address the key issues and improve the reliability of this 
sampling technique. 

The sensitivity of the abundance of clusters, both as a 
function of cluster mass and redshift, to the matter density 
has been used to estimate Om (see e.g., Bahcall & Fan, 
1998; Blanchard et al, 2000; Henry, 2000). However, the 
range in inferred values is broad, Om = 0.1 — 1, largely 
because of the exponentially important, but uncertain re- 
lation between cluster mass and x-ray temperature. For 
this reason, I do not include these measurements in my 
determination of the matter density. 

Finally, the value of the Hubble constant is important in 
converting physical densities to fractions of critical density; 
I adopt the value determined by the Hubble Key Project: 
h = 0.72 ± 0.07 (Freedman et al, 2001). The error is 
almost entirely due to systematics (the statistical error is 
only ah = ±0.02). The largest part of the systematic error 
budget is in the distance to the Large Magellanic Cloud. 
Possible dependence of the Cepheid period - luminosity 
relation on metallicity may also be an issue. 

3. THE CASE FOR O m = 0.33 ± 0.035 

Having adopted these physical measurements as input 
data, it is straightforward to deduce the fractions of critical 
density contributed by matter and by baryons, as well as 
the Hubble constant. But first, consider the consistency of 
these measurements. The physical baryon density {£l B h 2 ) 
is determined by BBN, CMB anisotropy and the power 
spectrum + Hubble constant: 

n B h 2 = 0.020 ±0.001 BBN 
= 0.022±°;°° 4 CMB 

= 0.022 ±0.011 Power Spectrum + H (6) 
These three, independent determinations of the baryon 
density are clearly consistent, giving one confidence in the 



case for a low baryon density (p B ~ 4 x 10~ 31 g cm -3 ). 
They involve very different physics - nuclear reactions 
when the Universe was seconds old, gravity-driven acous- 
tic oscillations when the Universe was around 400,000 yrs 
old, and the inhomogeneity in the distribution of matter 
in the Universe today - and thus also provide an impor- 
tant test the consistency of the big-bang framework and 
general relativity. 

Next, consider the ratio of the total matter density to 
the baryon density: 

Om/Ob = 7.2 ±2.1 CMB 

= 6.7 ± 3.1 Power Spectrum 

= 9.0 ±1.4 Clusters (SZ) + H 

= 8.7 ± 1.6 Clusters (x - ray) + H 

Again, all four measurements are clearly consistent and 
involve different physics - gravity driven acoustic oscilla- 
tions, inhomogeneity in the distribution of matter today, 
and cluster dynamics. 

Using standard techniques (and flat priors) I have con- 
structed the likelihood function for Om, Ob and h. From 
this, a posteriori probability distributions and credible 
ranges are calculated by marginalizing over one or two of 
these quantities in the usual way. The 1 — a ranges are: 

O s = 0.039± 0.0075 M = 0.33 ±0.035 

h = 0.69 ±0.06 (7) 

While the one-dimensional probability distributions are 
not perfectly Gaussian, the 68% and 95% credible ranges 
match pretty well with these la error flags. 

The value for the baryon fraction is largely driven by 
BBN and H Q ; alone they imply O s = 0.0385 ± 0.0077. 
Likewise, the Hubble constant is largely driven by its di- 
rect determination. The different scalings of the baryon- 
to-matter ratios with Hubble constant also provide impor- 
tant leverage, which can be seen in two ways: (1) the joint 
determination of h has a slightly smaller uncertainty than 
the direct measurement alone, ±0.06 vs. ±0.07; (2) if one 
arbitrarily doubles the uncertainty in the Hubble constant, 
h = 0.72 ± 0.14, and carries out the same analysis, the re- 
sults do not change dramatically: 

Ob = 0.040± 0.012 M = 0.33 ± 0.045 

h = 0.69 ±0.075 (8) 

To investigate the robustness of my estimates for the 
matter and baryon densities I have studied their sensitivity 
to the individual input data. To investigate sensitivity to 
scale-dependent bias on Hm/i I increased aa M h by a factor 
of 4. The central values for Om and Og were unchanged 
and their errors increased to ±0.044 and ±0.0076 respec- 
tively. One by one I doubled the errors for the other input 
data; the central values for Om and O B changed little and 
the errors grew to at most ±0.05 and ±0.008 respectively. 
One minor trend was observed: enlarging the error flags on 
both the X-ray and SZ cluster baryon fractions by a factor 
of 4, to significantly reduce the weight given to the clus- 
ter baryon fraction, decreased the estimate for the matter 
density by about la, Om = 0.29 ± 0.05. 

4. CONCLUDING REMARKS 
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Based upon present measurements of physical quanti- 
ties not tied to the distribution of light I conclude that: 
n M = 0.33 ± 0.035 and fl B = 0.039 ± 0.0075. While po- 
tential sources of systematic error remain, I have shown 
that these mass-density determinations are robust, and do 
not depend strongly upon any one measurement. 

The precision and reliability of these determinations 
should improve over the next decade. With Planck and 
MAP the uncertainty in both Q,Mh 2 and Vlsh 2 is likely 
to drop to the percent level (see e.g., Eisenstein et al, 
1999); the cluster and power spectrum derived quantities 
are likely to improve by at least a factor of two. The situa- 
tion with the Hubble constant is more difficult to predict. 

To illustrate the potential for improvement, I have re- 
duced the uncertainty in Q,Mh 2 and fish 2 to 1%, in Q,Mh 
to 5%, in Qb/^m to 20%, and in the cluster baryon frac- 
tions and Hubble constant by a factor of two. The pro- 
jected la uncertainties for VLm and SIb drop to ±0.012 and 
±0.0017 respectively. 

Already the physically based matter density has inter- 
esting implications. First, combining it with the CMB de- 
termination of the total mass/energy density, Oo = 1±0.05 
(Hanany et al, 2000; Netterfield et al, 2002; Pryke et al, 
2002), makes a very strong case for an additional compo- 
nent to the Universe (referred to as dark energy) that is 
smoothly distributed with Qx = 0.67 ± 0.06. This evi- 
dence for the dark energy bolsters significantly the direct 
evidence for dark energy from supernovae (Perlmutter et 
al, 1999; Riess et al, 1998). 



Second, the physically based matter density is about a 
factor of 2 larger than that determined from the mass-to- 
light ratios of clusters, Hm = 0.19 ± 0.04 (Carlberg et al, 
1997) and 0.16 ± 0.05 (Bahcall et al, 2000). This suggests 
that baryons in the cluster environment produce twice as 
much light as in the field. Going a step further, this im- 
plies the fraction of baryons that become stars is about a 
factor of two higher. The SDSS will determine mass-to- 
light ratio for a stacked sample of millions of field galaxies 
in five-color bands using weak lensing and will test this 
hypothesis by providing a reliable estimate of (M/L). 

Finally, in probing the nature of the dark energy, in- 
dependent knowledge of the matter density is crucial 
for breaking the degeneracy between the dark energy 
equation-of-state (wx = Px/px) and Q.m (see e.g., Maor 
et al, 2001; Wcllcr & Albrecht, 2001; Huterer & Turner, 
2001). The value deduced here is essentially independent 
of the nature of the dark energy (assuming a flat Universe), 
and the current one-sigma uncertainty would already help 
significantly to break the degeneracy in future determina- 
tions of wx ■ 
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